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1. Introduction 

The confinement of electromagnetic energy is an essential ingredient in studies of the 
quantum mechanical properties of light [U [21 [3] as well as in several applications 
involving the storage and manipulation of information [H ^J. In most circumstances, 
one requires light to be confined in small volumes over long periods of time. An optical 
cavity enables confinement of light through internal refiection on its boundaries [5j. 
The confinement is, however, only partial because some energy will always be lost to the 
surrounding environment. Therefore, the eigenmodes of these optical cavities — so-called 
quasi-normal modes — are characterized by a discrete set of complex eigenfrequencies, 
where the real part {u^) is proportional to the inverse of the wavelength of the confined 
light and the imaginary part (uj^^) is a measure of the temporal confinement of the 
wave inside the cavity. An important figure of merit is the quality factor Q of these 
modes, which is usually defined as the temporal confinement of the energy normalized 
to the frequency of oscillation, such that represents the fraction of energy lost 
in a single optical cycle. The quality factor can be calculated as Q = uj'/{2uj"). 
Using specific fabrication techniques, dielectric resonators that exhibit extremely high 
quality factors have been realized. Experimentally, quality factors up to 8 x 10^ 
have been measured in dielectric microsphere resonators [6j and larger than 10^ 
in dielectric toroid microcavities on a chip [7J. Cavities with high quality factors 
in combination with small mode volumes are extremely interesting for applications 
involving cavity quantum electrodynamics [HI [9j, such as the recent developments 
to integrate optical microresonators into atom chips [lOj. In these applications, 
it is important to have a small vacuum region in which atoms interact with the 
electromagnetic modes. Unfortunately, dielectric cavities are fundamentally limited in 
size because it is impossible to efficiently store light in volumes with dimensions smaller 
than the wavelength of the confined mode |im [12]. One attempt to overcome this 
limitation and to miniaturize the mode volume of the confined light is the development 
of surface plasmon polariton cavities [13j. Here, however, the temporal confinement is 
severely reduced by dissipation in the metals. 

In this paper, we show that — using the formalism of transformation optics — 
alternative designs of subwavelength optical cavities exist. This approach is based on the 
equivalence between Maxwell's equations in vacuum, expressed in a curved coordinate 
system, and Maxwell's equations inside a nontrivial material with specific permittivity e 
and permeability ji p!4l [T5l [T6] . Although this mathematical equivalence was known for 
quite some time [TTl [HI [IHl [20] , it was only recently proposed to effectively realize such 
coordinate transformations with the use of metamaterials [I4[[15]. Transformation optics 
has demonstrated its huge potential through various proposals of novel optical devices 
that manipulate the electromagnetic beams in unconventional ways [21j. Amongst 
many others, transformation optics has already been used to design perfect lenses 
p2l [231 [24l [25] . beam and polarization manipulators j26l [27], super scatterers [28j, 
invisibility cloaks [HI [151 [30l |31j , and devices implementing other optical illusions 
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[321 [331 134] . Moreover, the general four-dimensional formulation of transformation optics 
[35] allows for applications that also involve the time coordinate, such as a frequency 
convertor [36l |37l [38], a laser pulse analogue of Hawking radiation [39l [40l HTJ, an 
electromagnetic analogue of Schwarzschild-(anti-)de Sitter spacetime [12], or a spacetime 
cloak jSj. 

The specific permittivity and permeability tensors required to implement devices 
designed using the techniques of transformation optics usually do not exist in nature 
and must, therefore, be achieved with the aid of metamaterials. These man- 
made artificial materials derive their electromagnetic properties from subwavelength, 
appropriately designed constituents In particular, metamaterials can be made 

with a negative pemittivity and permeability at the same frequency. These so-called left- 
handed materials exhibit peculiar phenomena such as negative phase velocity, negative 
refraction, and inversed Doppler effect [45j. Thanks to the ability to compensate the 
phase of electromagnetic waves inside left-handed materials, they generate the possibility 
of perfect imaging [46j and miniaturization of optical devices beyond the diflFraction limit 
[171 UHl [19]. More recently, considerable interest has been devoted to metamaterials 
with other electromagnetic functionalities such as artifical chirality [50l [5lj , a classical 
analogue of electromagnetically induced transparency [52[ [53l l54] , or the enhancement 
of quantum phenomena [55] . 

In this paper, we will investigate several approaches to design an optical cavity 
within the framework of transformation optics. In section [2} we will derive the dispersion 
relation of a cavity based on a radial coordinate transformation and we will apply 
this dispersion relation to calculate the bound modes of a cavity in which the radial 
coordinate is transformed under a hyperbolic function. In section [3j we introduce a 
folded coordinate transformation, which results in a left-handed cavity characterized 
by a continuum of eigenmodes with an infinite quality factor. Finally, in section [4| we 
introduce a combined transformation on the radial and azimuthal coordinate, which is 
used in section [5] to show how perfect confinement can be achieved inside a cavity made 
of right-handed materials. 

2. A hyperbolic cavity 

2.1. The dispersion relation of a cylindrical cavity based on a radial coordinate 
transformation 

We start by calculating the bound modes of the system shown in figure [T} which 
consists of a hollow cylinder bound by the radii p = Ri and p = R2. To obtain 
the eigenmodes of this system, we solve Maxwell's equations inside each region and 
match the solutions using the proper boundary conditions. Without loss of generality, 
we will consider the time-harmonic solutions that are polarized along the z-axis (TE- 
polarization) : E(r, t) = E[p^ cp) exp(— icjt) 1^. 

Inside the empty regions. Maxwell's equations combine into the free-space Helm- 
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Figure 1. An infinite hollow cylinder (II) with inner radius Ri and outer radius 
surrounded by vacuum (III), includes an inner region (I) where electromagnetic fields 
might be trapped. The medium that performs the coordinate transformation — called 
"transformation-optical medium" — is situated in region (II) 



holtz equation for the electric field. The angular variation of the electric field is a 
sum of imaginary exponentials, characterized by the indices ui and whereas the 
radial dependence satisfies the cylindrical Bessel equation. To simplify the boundary 
conditions in the next step, we will use the Bessel functions and Y^y in region (I) and 
the Hankel functions hI^"^ and hI^"^ in the surrounding region (III), because the latter 
can be interpreted as incoming and outgoing cylindrical solutions of the Bessel equation: 

Ei{p, 0) = [A^j Juiikop) + Y^^{kop)] exp(iz/i0), (1) 

Euiip. ct>) = [a„, Hillikop) + i^.,, Hl^likop)] exp(iz.ni0), (2) 

where Aj,^^ S^j, C^^jjj, and Dj,^^^ are complex integration constants and ko = uj/c 
represents the vacuum wavenumber. 

In order to calculate the solutions in the transformation-optical region (II), we 
need to insert the constitutive parameters for the material in Maxwell's equations. In 
transformation optics, these constitutive equations can be derived from the coordinate 
transformations required to impose a specific pathway onto the electromagnetic fields 
[lllIlSlII^. In this section, we consider the case of an arbitrary radial transformation, 
leaving the azimuthal angle and the z-axis unchanged: (p, 0, z) is transformed into 
{p\ (\)\ z') such that 

p' = /(p), (3) 
<t>' = 4>, (4) 
z' = z. (5) 

This distortion of the radial coordinate in vacuum has the same effect on the 
electromagnetic radiation as if it were propagating in a medium with the following 
nonzero components of the permittivity and permeability tensors: 
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_ 6 _Pfip). (Q) 

. ^ . ^ /(p)r(p) 

where prime denotes differentiation. 

Using the constitutive equations = ji^ji^-H^ and = e^e^-E^^ we can insert 
these parameters in MaxweU's equations and combine them into the foUowing equation 
for the electric field in region (II): 



(7) 



The solutions of equation ([T]) will of course keep their harmonic azimuthal character. 
As to the radial part of this equation, it can be reduced to the same Bessel equation as 
in region (I) and (III), but in the variable p' = /(p). As a result, the solutions inside 
the transformation-optical region (II) are given by 

En{p. (/>) = J^Mip)) + Y^Mip))] exp(iz.„0), (8) 

where, once again, Fj^^^ and Gj^^^ are complex integration constants. 

This solution can now be matched to the solutions in (I) and (III) [equations ([T]) 
and ([2])], using the appropriate boundary conditions. Obviously, both the electric and 
magnetic fields should be periodic in 0, or E{p, 0) = E{p, 27r)|fl This condition is fulfilled 
when Ui = rrii^ with rrii G Z for all regions i. Furthermore, we consider only those modes 
whose amplitude is finite, which implies that we should reject the Bessel function in 
region (I), because it has a singularity at the origin. In region (III), on the other hand, 
we impose Sommerfeld's radiation condition, expressing that no energy is fiowing in from 
infinity. We should therefore drop the second Hankel function i^m , which represents 
such an incoming wave. 

The dispersion relation is now found by imposing the continuity of the tangential 
components of the electric (E^) and magnetic fields (H^) at the boundaries between 
the regions (I), (II), and (III). Firstly, since the boundaries p = Ri and p = R2 do 
not depend on 0, the azimuthal mode numbers rrii must be the same in each region. 
Secondly, we find the following set of four independent equations, in which we already 
eliminated the angular parts: 

JUkoRi) = JUkofiRi)) + Y^ikofiRi)), (9) 

.UhR,) = l^j;^{k,f{R,)) + Gm ^^y;(fco/(i?i)), (10) 

JrCi JrCi 
Fm Jm{hfiR2)) + Gm ^^(fco/li^s)) = H^\koR2), (11) 

Fm ^J'Jkof{R2)) + Gm Z^y;(fco/(i?2)) = Cm H'i'\koR2), (12) 

where ^4^, F^, G^, and are complex integration constants. The equations involving 
the magnetic field were simplified using the relation f\p) / = f{p)/ p, which is derived 

:|: This also implies that the magnetic field is periodic in this direction. 
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from equation ([6]). Setting the determinant of this set equal to zero generates the 
dispersion relation of the system and determines the eigenmodes of the cavity. Note 
that this relation is valid for any cavity of the type shown in figure [l} implementing a 
radial coordinate transformation p' = f{p) between Ri and i?2. 



2.2. The confined modes inside a hyperbolic cavity 



A cavity ideally confines the electromagnetic energy in a small (subwavelength) region 
of space for a very long time. In terms of electromagnetic and physical space, this 
could be achieved by mapping some large domain of the electromagnetic space onto a 
much smaller region in physical space. Such a transformation can be constructed with 
a hyperbolic function that grows to infinity in a finite point. We will therefore consider 
a device as shown in figure [T| in which the following radial coordinate transformation is 
implemented between Ri and R2: f : i?2] [Ri^ oc] : p p', where 

f{Ri) = Ru (13) 
f{R2) = 00. (14) 
The matching at the inner boundary enables a smooth transition of the waves. 



Since there cannot be anything "beyond infinity," equation (14) ensures that the 



electromagnetic energy cannot escape this device. Such a coordinate transformation 
is illustrated in figure [2](a), where we consider the transformation given by /(p) = 
Ri{Ri — R2) I (p— ^2)- The Cartesian coordinate lines in physical space become denser as 
we approach the outer radius R2. A similar transformation has been proposed to design 
a matching layer in order to improve the efficiency of numerical software algorithms 
[56] . The values of the material parameters required to implement this transformation 
in physical space are shown in figure [2](b). 

The determinant of equations ([9])-(12) can be calculated in the limit of /(i?2) 
oc, using equation (13). We find that the hyperbolic map does not confine any 




0.6 0.7 0.8 0.9 1.0 



Figure 2. (a) The underlying coordinate transformation of a hyperbolic cavity, defined 
by /(p) = Ri{Ri — R2)/{p — R2)' The region i?2] in the physical radial coordinate 
p covers the region [i?i,+oo] in the electromagnetic radial coordinate p' . (b) The 
material parameters that implement this hyperbolic transformation. In the limiting 
case of f{R2) 00 the radial components of e and p become zero, while the other 
parameters grow to infinity at the outer boundary. 
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Figure 3. The confined modes with azimuthal mode number m = 3 of a perturbed 
hyperbolic map, in which we cut out a rim AR of the outer boundary i?2, such that 
the outer value Ro^t = R2 — Ai? is not mapped onto infinity, but instead takes the 
value f{R2 — AR) = 20. (a) A contour plot in de complex frequency plane indicating 
the solutions of the dispersion relation, (b) The two-dimensional plot of the electric 
field corresponding to the first solution of the dispersion relation. 

electromagnetic modes. In contrast to what is mentioned in reference ^57j, we find 
that, independent of the azimuthal mode number m, the dispersion relation can only be 
satisfied if ko = 0, i.e., the static solution. This result fits with the intuitive idea that in 
this configuration an electromagnetic wave travels an infinitely long time to reach the 
outer boundary of the cavity. Therefore, no standing wave can be created in the cavity: 
its structure does not permit a reflected wave at p = R2. 

The previous physical interpretation implies that a perturbed version of the 
hyperbolic map — in which the material parameters do not grow to infinity — should 
exhibit confined modes. This is indeed confirmed by the numerical evaluation of the 
dispersion relation, whose solutions are shown in figure jS) The number of solutions 
increases as AR decreases, or equivalently as /(i?2 — AR) approaches oc. The quality 
factor Q increases at the same pace. As shown in figure [3| the cavity enables 
subwavelength confinement of electromagnetic energy; the first solution, for example, 
lies at /coi?2 = 0.51 — 7.0 x 10~^ i, which corresponds to a free space wavelength Aq that 
is more than 10 times larger than the outer radius of the cavity R2. The quality factor 
of this mode is Q = 3.6 x 10^. 

In figure [3](b), we show a two-dimensional plot of this mode inside the 
cavity. The reader should note how the field is almost completely located in the 
transformation-optical medium, which sounds reasonable since this medium contains the 
electromagnetic interval [i?i,+oc], whereas the inner disk (region I) only occupies the 
electromagnetic interval [0, inside the transformation-optical region the wavelength 
of the electric field becomes smaller towards the outer radius R2 due to the increasing 
values of the material parameters inside the medium. 

Judging from these results, the imperfect hyperbolic design might seem to be the 
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subwavelength optical cavity we are looking for, having good confinement in arbitrarily 
small dimensions. We should, however, look back at the materials with which it is 
implemented in figure [2](b). Firstly, we note that a traditional whispering gallery 
cavity, made entirely from these high-index materials, also has subwavelength modes. 
Secondly, the wavelength is becoming extremely small within the device so that it 
is practically impossible to use the mean-field approximation when determining the 
material parameters. In the subsequent sections, we will present subwavelength cavities 
in which this is no longer the case. 



3. The perfect cavity 



In the example of the hyperbolic cavity, the entire electromagnetic space was mapped 
onto a finite region of physical space. We can, however, consider a cavity from a cloaking 
perspective and design a device that cloaks away the volume surrounding the device, 
instead of the volume inside the device [58j. Such a device should smoothly guide the 
electromagnetic waves in the cavity so that they never penetrate the outer boundary. 
The eflFect of such a transformation is shown in figure [4](a). Since we want to cloak 
away region (III), we will use a radial coordinate transformation that maps the physical 
coordinates (p, 0, z) onto the electromagnetic coordinates (p^, 0', z'). To achieve perfect 
cloaking of region (III) from the viewpoint of region (I), the radial transformation 
function has to satisfy the following boundary requirements: 

= (15) 

/(i?2) = 0. (16) 

As before, the actual shape of the function has no implications on the cavity's 
performance. Transformation functions that satisfy these boundary conditions have also 
been studied in combination with traditional invisibility cloaks, giving rise to so-called 
anti-cloaks [59l [60j . 

The modes of the present cavity are the solutions of equations ([9])-(12), where we 
now have to insert f{Ri) = Ri and /(i?2) = 0. These equations now become 

Am Jm{koRi) = Fm Jm{koRi) + Gm Ym{k^Ri) , (17) 

Am JLikoRl) = Fm JmiW + Gm Y^^ikoRl), (18) 

Fm lim Jm{kox) + Gm lim Ym{kox) = Cm H^^\koR2), (19) 



Frr 



lim 

x^O 



X 

R2 



Jmikox) 



x^O 



lim 



X 



= CmHi^\koR2). (20) 



These limits should be handled with care, since they contain indefinite expressions like 
X oc. Assuming the azimuthal mode number m 7^ 0, these limits can be unambiguously 
evaluated: 



lim Jm{kox) 0, 



x^O 



limYm(kox) = — oc. 



(21) 
(22) 
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Figure 4. (a) The coordinate lines that are generated by a radial coordinate 
transformation implementing a perfect cavity, defined by f{p) = -^^r^{p — R2)' 

(b) The material parameters required to materialize this coordinate transformation. 

(c) The electric field mode profile of a confined mode inside the cavity. The wavelength 
of this mode is much larger (factor 100) than the outer radius of the cavity. 



lim 



lim 



X 



^2 



= 0, 



X 



-^Y^ikox) 



+ OC. 



(23) 



(24) 



We can now reinsert these limits in equations (19)-(20) and we find that this set only 
has solutions if Gm = and Cm = for all azimuthal mode numbers m 7^ 0, whereas 
there are no requirements on F^. Taking this into account, equations (17) and (18) 
become 

Am JmikoRl) = Fm Jm{hRl), (25) 
Am J'mihRl) = Fm JmWl). (26) 

hence Am = Fm- This set imposes no constraints on /cq, which means that the cavity 
supports a continuous spectrum of modes, even if the wavelength is larger than the 
characteristic dimensions of the cavity. These modes are perfectly confined, since D is 
equal to zero: there is no radiation escaping into region (III). The quality factor Q is 
infinite and, as a consequence, the complex part of the frequency (cj") should be zero. 

We are now able to plot the solutions of the perfect cavity. We can choose any real 
free-space wave vector ko and plot the solutions, using equation ([s]). In figure [4](c), we 
plot a mode for which koRi = 0.01. The field's variation inside the cavity depends on 
the chosen transformation function /(p). One can make well-considered choices for this 
function / to enhance the field distribution inside the transformation-optical medium. 

We observe a completely different mechanism of confinement as compared to the 
hyperbolic map. Generally, a wave can be confined inside a cavity if one round trip 
(approximately the cavity's circumference) equals an integer number / of the mode's 
wavelength inside the cavity: 27ra ^ IX [llj. The perturbed version of the hyperbolic 
map reduces the wavelength of an electromagnetic mode to a very small number at the 
outer boundary, thus fulfilling the condition. In the perfect cavity, however, the phase 
shift vanishes completely and / = 0. 
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The reduction of equations (17)- (20) to the trivial equations (25)- (26) was only 
possible when we assumed the azimuthal mode number m 0. A mode without 
azimuthal momentum cannot be confined within this cavity. Physically, this can be 
understood since such a mode has a purely radial wave vector and in the absence 
of azimuthal propagation it cannot be deflected to the left or to the right inside the 
transformation-optical region. 

The material losses are high due to the fact that the transformation-optical medium 
is made of left-handed materials, as shown in flgure |4](b). Although the material 
parameters strongly depend on the choice of the transformation function /(p) — through 
equations ^ — one can prove that any transformation-optical medium that satisfles 



equations (15)-(16) will have a region in which all components of the permitivity and 
the permeability tensors are negative. This is analogous to the perfect lens [46l EI] — 
another example of a folded map [16j — which also requires a left-handed response. 

In the last sections of this article, we will derive a method to overcome this limitation 
and demonstrate how it is possible to design a cavity with right-handed material 
parameters only. But let us flrst introduce the idea of cavities based on azimuthal 
coordinate transformations. 

4. Azimuthal coordinate transformations 

^.7. The dispersion relation in case of an azimuthal coordinate transformation 

In this section, we investigate transformation-optical cavities, as shown flgure [T} in 
which the transformation also involves the azimuthal coordinate (j). We will consider a 
transformation deflned by 

p' = f{p), (27) 
4>' = gi4>), (28) 
z' = z. (29) 

and, once again, look at solutions of Helmholtz' equation with linear polarization along 
the ;^-axis. It can be shown that such a transformation can be implemented with 
materials whose components are 

m 



f{p)f'{p) u,-. 
^ z=P z = 9 {4>). 

in which f\p) denotes differentiation of /(p) with respect to p and 5''(0) denotes 
differentiation of g{(j)) with respect to (j). The wave equation of such a medium is 
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whose solutions are given by 

E,,{p, 4>) = [F, J^Mip)) + Y^Mip))] exp{iuug{(l>)). (32) 

Here again, the cylindrical symmetry leads to the quantization of the azimuthal mode 
number uu rriu'. 

""(*^' = 9(2.) -9(0)' 

with k G Z. However, unlike the dispersion relation derived in previous sections 
where the mode numbers in the different regions were identical, a general azimuthal 
transformation will scramble the azimuthal momenta of the solutions in the different 
regions. One single azimuthal mode exp{imiig {(/))) in the transformation-optical region 
(H) will excite multiple modes in the vacuum region, and vice versa: 

exp(imii^(0)) = ^Cmi exp(imi0), (34) 

mi 

where the coefficients Cmi are given by 

1 f^"" 

Cmi = ^ / exp(i (mii^(0) - mi0))d0. (35) 
2vr Jo 



The Fourier series expansion in equations (34)-(35) is possible since exp^imug {(/))) is 
a periodic function of the azimuthal coordinate (j) with period 27r, as can be seen by 
substituting equation (33) into equation (32). In the surrounding vacuum region (HI), 
the same condition on the azimuthal coordinate applies, Cmi = C'^nr 

We will calculate now the dispersion relation and restrict the analysis to linear 
azimuthal transformations g{(l)) = acj)^ where a is a real number. Using equation (35), it 
can be shown that in this case a single azimuthal mode number mi = mm = mi in the 
vacuum regions will match with a single mode number mu = m2 in the transformation- 
optical region, where these mode numbers are related by 

1712 = — . (36) 
a 

In general, the angular mode number m2 will not be an integral number. The dispersion 
relation is then similar to the one that corresponds to a single radial coordinate 
transformation ([9])-(12) and is generated by the following set of equations: 

Am, JmAkoRl) = JmM{Rl)) + ^m, Y^^ih f {Rl)) , (37) 

Am, J'm,{k,Ri) = Fm, {Rl)) 



+Gm. ^V^K..(fco/(i?i)), (38) 



f{Ri)a^ 
Ri 

Fm, JmM{R2)) + Gm, {R2)) = ^m, ^S(fcoi?2), (39) 



Fm. ^-^JUk,f{R2))+Gm. ^^^K..(W(i?2)) 

= Cm, HS{koR2)- (40) 
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Figure 5. (a) Contour plot of the dispersion relation of the confined modes whose 
azimuthal mode number in region (I) equals mi = 5. The cavity is defined by equation 
( [42| ). (b) A density plot of the electric field distribution inside the cavity, corresponding 
to the fourth solution in (a), at ko = 13 — 5.7 x 10~^ i. 



The additional factors a = g^c/)) in equations (38) and (40) originate from the fact that 



m _ r(p)f(p)^-^^) ^ m (41) 



Pf'ip) 



4-2. Cavities based on a single azimuthal transformation 

It is instructive to have a look at the modes of a cavity defined by an exclusively 
azimuthal transformation, for instance, 

g{<P) = 150, (42) 

fir) = r. (43) 

As shown in figure [5](a), there are several confined modes within these cavities. A 
corresponding mode profile is shown in figure [5](b). Very much in correspondence with 
the hyperbolic map on the radial coordinate in section [2} we find that the quality factor 
of the solutions increases as the optical path length inside the cavity increases. In 
contrast to the hyperbolic map, however, the number of subwavelength solutions does 
not drastically increase as we increase the optical path length. 

Another intriguing example is the collapse of the azimuthal coordinate, i.e., all 
angles are transformed on one and the same angle {g{(l)) = (/>o), inside a full cylinder (no 
vacuum region). Although the corresponding material parameters are extremely exotic 
(zero and infinity) and thus not useful for practical applications, this setup is interesting 
for theoretical reasons. Obviously, there will be no quantization of the azimuthal mode 
number m, and the continuity relations imply that 

/ + CX) 
F{m) JmikoR)dm = C H'iXhR), (44) 
-oo 

/ + CX) 
F{m) J^{koR)dm = C H^^\koR). (45) 
-OO 
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which immediately translates into 

/ + CX) 
F{m) J^{koR)dm = 0. (46) 
-OO 

This equation has a solution for any /cq. The cavity thus confines light at every 
wavelength. These two examples clearly show the difference between a radial and 
an azimuthal transformation. The former changes the radial coordinate p, which 
automatically alters the quantization of /cq, whereas the latter manipulates 0, which 
changes the azimuthal mode number m and thus only indirectly infiuences the 
quantization of ko through the dispersion relation. 



5. Perfect confinement in a right-handed cavity 



Azimuthal transformation optics can be very valuable when used in combination with 
a nontrivial radial transformation. To demonstrate this, we show here how the addition 
of an azimuthal transformation can be used to generate a cavity with right-handed 
material parameters in which there are no radiation losses. Let us consider the radial 
transformation of a perfect cavity, in combination with an azimuthal transformation 
that inverts 6: 



fift) 

9i<l>) 



Ri 



i?2 



2) 



(47) 
(48) 



implemented between the radii p = Ri and p = R2. 

The material parameters will be the same as those of a perfect cavity. However, 
due to the inversion of (j) an additional sign reversal will make all material parameters 
positive. The dispersion relation of this cavity is then given by the following set of 
equations: 



Am Jm{koRi 

R 

i-xl 
R2 



m lim J_rn{koX) + G 



F-m lim 

x->0 



J-m{koRl) + G-m y-m{koRl) , 

{-l)J'_^{koRi) + G.^ {-l)YLJkoRi), 
{kox) = C„^H^\koR2), 
i-x). 



x-)-0 



G-m lim 



x->0 



R2 



-YL^{kox) 



G^H'^'\koR2). 



(49) 
(50) 
(51) 



(52) 



Following the same argumentation as for equations rtl9h and (20), we find that equations 



(50) and (52) can be solved if G_m = and Cm = for all azimuthal mode numbers 
m 7^ 0, without any restrictions on the values of and ko. We can reinsert this in 
the boundary conditions at p = i?i : 

Am Jm{koRl) = F_m J-m{koRl), (53) 



Am JLikoRl) = F_m {-l)f_m{koRl)- 



(54) 
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Figure 6. (a) The grid lines of a right-handed cavity in which subwavelength modes 
can be confined in the absence of radiation losses. These grid lines correspond 



to a transformation defined by equations (61)- (63) (b) The material parameters 
required to materialize this coordinate transformation. These material parameters 
were considerably simplified by choosing Ri = O.7R2 (c) The electric field mode profile 
of a perfectly confined mode inside the cavity (m = 1, Aq ~ 3AR2). 



Using the identity = (— l)^J^(x), it is clear that the set can be solved for all 

frequencies for which Jm{koRi) = or J^(fcoi?i) = 0. Depending on the angular mode 
number m, A^^ then equals or — 

The eigenfrequencies of this cavity are defined by the zeros of the Bessel's function 
or its derivative: 

ko=jm,n/Rl^ (55) 

ko=fmjRi^ (56) 

where jm,n and ,^ are the nth solution of Jm{x) = and J^{x) = respectively. The 
eigenfrequencies, therefore, cannot be chosen at will. To overcome this limitation we 
modify the design by replacing the vacuum in the inner region (I) by a transformation- 
optical material that maps the radial coordinate onto a larger one, i.e., f{Ri) = fn^ > 
Ri. Inside the inner region we implement the transformation given by /(p) = fn^p/Ri- 
Obviously, the transformation in region (II) should also map Ri onto /r^. This 
modification allows us to design the cavity to have perfectly confined modes at arbitrary 
frequencies, since the resulting dispersion relation is given by: 

Jm{k{)fR^) = F-m J-m{k{)fR^), (57) 

Am J'Mr,) = F_m {-1)J'_Mr,). (58) 
This cavity has solutions for fco/i^i = or fco/i^i = jm.n- Equivalently, we can write 

R2 R2 jm,n 



i?2 R2 Jm,n 



(59) 
(60) 



In figure loFa), we plot the grid lines of such a cavity in which the inner radius is 



15 



mapped onto a larger value: Jr-^ = i?2- The underlying coordinate transformations are 

R2 

Ri 



m = ^P, (61) 



in region (I), and 



R2 



f{p) = > \ - Rl (62) 

^ ^ ^yRl -Rr 

gi4>) = -0, (63) 

in region (II). The thick red hne clearly indicates that the coordinates in the inner region 
are not continuously guided into region (II). The necessary condition for reflectionless 
transformation media is not valid since at the interface p = Ri the coordinates of region 
(II) cannot be matched with those of region (I) through a combination of rotation 
and displacement [62 J. This cavity, therefore, only confines light at discrete resonance 
frequencies. This gives a geometrical explanation of the discreteness of the solutions as 
given by equations (59)- (60). Figure [6](b) shows the material parameters of this cavity. 
Three elements of the material tensors are simplified to the vacuum values thanks to 
the particular choice of Ri = 0.7i?2. Finally, in figure [6](c), we plot the electric field of 
a perfectly confined mode inside this cavity. There are no fields outside the cavity and 
the cavity is subwavelength (Aq = 3.4i?2). 



6. Conclusion 



In this paper, we have discussed designs of electromagnetic cavities based on 
transformation optics. We derived the dispersion relations of cavity structures based 
on radial and azimuthal coordinate transformations and applied those to calculate 
their bound modes. Some of these transformations enlarge the optical path length 
inside the cavity, whereas others are based on a folding of the electromagnetic space. 
Finally, we have shown how the combination of radial and azimuthal transformations can 
eliminate the left-handedness of the perfect cavity, while preserving its most important 
characteristics: confinement of electromagnetic modes with unlimited quality factor due 
to radiation losses, even if the wavelength is larger than the dimensions of the cavity. 
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